o 



(N 



Second- Order Theory for Iteration Stable 

Tessellations 

by Tomasz Schreiber and Christoph Thale 



^ ! from Toruri and Fribourg 



Abstract 



This paper deals with iteration stable (STIT) tessellations, and, more generally, 
Ph . with a certain class of tessellations that are infinitely divisible with respect to itera- 

tion. They form a new, rich and flexible class of spatio-temporal models considered 
in stochastic geometry. The martingale tools developed in [17] are used to study 
second-order properties of STIT tessellations. Firstly, a general formula for the vari- 
ance of the total surface area of cell boundaries inside a convex observation window 
is shown. This general expression is combined with tools from integral geometry 
to derive explicit exact and asymptotic second-order formulae in the stationary and 
isotropic set-up, where a family of chord-power integrals plays an important role. 
■ Also a general formula for the pair-correlation function of the surface measure is 

OH ■ found. 
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h : 1 Introduction 

Iteration stable random tessellations (or mosaics), called STIT tessellations for short, form 
a new model for random tessellations of the rf-dimensional Euclidean space and were for- 
mally introduced in [10, 11, 12, 13]. They have quickly attracted considerable interest in 
stochastic geometry, because of their flexibility and analytical tractability. They clearly 
show the potential to become a new mathematical reference model beside hyperplane and 
Voronoi tessellations studied in classical stochastic geometry. Whereas much research in the 
last decades was devoted to mean values and mean value relations, modern stochastic ge- 
ometry focusses on second-order theory and distributional results, see [1, 2, 3, 4, 5, 6, 7, 14] 
to mention just a few. 

To introduce the non-specialized reader to the subject, we briefly recall the basic construc- 
tion of STIT tessellations within compact convex windows W C R d with interior points. 
To this end, let us fix a (in some sense non-degenerate) translation-invariant measure A 
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on the space of hyperplanes. Further, let t > be fixed and assign to the window W a 
random lifetime. Upon expiry of its lifetime, the primordial cell W dies and splits into two 
sub-cells separated by a hyperplane hitting W, which is chosen according to the normalized 
distribution A. The resulting new cells are again assigned independent random lifetimes 
and the entire construction continues recursively until the deterministic time threshold t 
is reached (see Figure 1 for an illustration). In order to ensure the Markov property of the 
above construction in the continuous-time parameter t, we assume from now on that the 
lifetimes are exponentially distributed. Moreover, we assume that the parameter of the 
exponentially distributed lifetimes of individual cells [c] equals A([c]), where [c] stands for 
the collection of hyperplanes hitting c. In this special situation, the random tessellation 
constructed by the described dynamics fulfils a stochastic stability property under the op- 
eration of iteration of tessellations, and whence is indeed a STIT tessellations. We refer 
to Section 2 below for more details. 

In [17] we have introduced a new technique relying on martingale theory for studying these 
tessellations. One feature of this new approach is that it allows to investigate second-order 
parameters (i.e. variances) of the tessellation, which were out of reach so far and are in 
the focus of the present work. Based on a specialization of our martingale technique, we 
calculate in Section 4 the variance of a general face-functional and as a special case we find 
the variance of the total surface area of cell boundaries in a bounded convex window. The 
resulting integral expression can be explicitly evaluated in the stationary and isotropic 
case by applying an integral-geometric transformation formula of Blaschke-Petkantschin 
type, which is also developed in this paper, see Section 5. For the particular case of space 
dimension 3, an exact formula without further integrals is found. Another important task 
in our context is to determine for fixed terminal times t the large scale asymptotics of 
the afore calculated exact variance for a family of growing compact and convex windows 
with positive volume. Relying again on techniques from integral geometry, we will be able 
to determine asymptotic variance expressions, leading - most interestingly - in dimension 
d — 2 to a result of very different qualitative nature compared with space dimensions 
> 3, where certain chord-power integrals, known from convex and integral geometry, will 
reflect the influence of the geometry of the observation window. In dimension d = 2 we 
will see that, in contrast to the described situation for higher space dimensions, the shape 
of the window does not play any role and only its area enters our formulas, see Section 
6. We also derive an explicit expression for the so-called pair-correlation function of the 
random surface measure for arbitrary space dimensions, see Section 7, generalizing thereby 
recent findings from [22], which are based on completely different methods. This function 
is a a commonly used tool in spatial statistics and stochastic geometry to describe the 
second-order structure of a random set and it describes the expected surface density of the 
tessellation at a given distance from a typical point. 

We would like to point out that the second-order theory developed in this paper is funda- 
mental for our further work on STIT tessellations [18, 19, 20] and that its extended version 
[16] is available online. 

In this paper we will make use of the following notation: 
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Figure 1: Realizations of a planar and a spatial stationary and isotropic STIT tessellation 
(kindly provided by Joachim Ohser and Claudia Redenbach) 



• Br = Bji(p) is the d- dimensional ball around the origin with radius R > 0. 

• Kj := Vol j(B{) is the volume of the j-dimensional unit ball, jKj its surface area. 

• The uniform distribution on the unit sphere Sd-i in M. d (normalized spherical surface 
measure) is denoted by Vd-x- 



2 Construction and properties of the tessellations 

Let A be a non-atomic and locally finite measure on the space H of hyperplanes in the 
d- dimensional Euclidean space ~R d . Further, let t > and W C M. d be a compact convex 
window with interior points in which our construction of a random tessellation Y(tA, W) 
is carried out. In a first step, we assign to the window W an exponentially distributed 
random lifetime with parameter A([W]) where [W] := {H G "H, H D W ^ 0} stands 
for the collection of hyperplanes hitting W. Upon expiry of its lifetime, the cell W dies 
and splits into two polyhedral sub-cells W + and W~ separated by a hyperplane in [W], 
which is chosen according to the law A(-)/A([W]). The resulting new cells W + and W~ 
are again assigned independent exponential lifetimes with respective parameters A([W + ]) 
and A([W~]) (whence smaller cells live stochastically longer) and the entire construction 
continues recursively until the deterministic time threshold t is reached (for an illustration 
see Figure 1). The cell-separating (d — l)-dimensional facets (the word facet stands for 
a (d — l)-dimensional face here and throughout) arising in subsequent splits are usually 
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referred to as (d — l)-dimensional maximal polytopes (or I-segments for d — 2 as assuming 
shapes similar to the letter T). 

The described process of recursive cell divisions is called the MNW-construction in honour 
of its inventors in the sequel and the resulting random tessellation created inside W is 
denoted by Y(tA, W) as mentioned above. The random tessellation Y(tA, W) has the 
following properties (see [13] for detailed proofs): 

• Y(tA, W) is consistent in that Y(tA, W) n V = Y(tA, V) for convex V C W and thus 
Y(tA, W) can be extended to random tessellation Y(tA) on the whole space M. d . 

• If A is translation- invariant, Y(tA) is stationary, i.e. stochastically translation in- 
variant. If, moreover, A is the unit-density isometry- invariant hyperplane measure 
A iso , then Y(tA iso ) is even isotropic, i.e. stochastically invariant under rotations wrt. 
the origin. 

• Y(tA, W) is iteration infinitely divisible with respect to the operation EH of iteration 
if tessellations for any compact convex W C M d . This is to say 

Y(tA, W) = m(Y((t/m)A, W) EE! • ■ ■ EE! Y((t/m)A, W)), m = 2, 3, . . . , 

explained in detail in [17]. Because of this property we call Y(tA, W) an iteration 
infinitely divisible MNW-tessellation. In addition, if A is translation-invariant, Y(tA) 
is stable under the operation EB, which is to say 

Y(tA) = m(Y((t/m)A) EB • • • EB Y((t/m)A)), m = 2, 3, . . . . 
For this reason, Y(tA) is called a random STIT tessellation in this case. 

• In the stationary set-up, the surface density, i.e. the mean surface area of cell bound- 
aries of Y(tA) per unit volume equals t. 

• STIT tessellations have the following scaling property: 

tY(tA) = Y{A), 

i.e. the tessellation Y(tA) of surface intensity t upon rescaling by factor t has the 
same distribution as Y(A), the STIT tessellation with surface intensity 1. 

3 Background material 

In this section we recall a few facts from [17], wich are going to be crucial for our argu- 
ments below. Firstly, it follows directly from the MNW-construction of Y(tA, W) that, 
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in the continuous-time parameter t, this is a pure jump Markov process on the space of 
tessellations of W , whose generator is L := h A . w with 

LF(Y)= [ [F(YU{f})-F(Y)]A(dH) (1) 

J\W] ^„_T7„„rr, 



/eCeiis(ynff) 

for all F bounded and measurable on space of tessellations of W. Similar to the approch 
taken in [17], the general theory of Markov processes can now be used to construct a class 
of martingales associated with iteration infinitely divisible MNW-tessellations or, more 
specifically, STIT tessellations. Indeed, for bounded measurable G = G(Y,t), considering 
the time-augmented Markov process (Y(tA,W),t) t > and applying standard theory, see 
Lemma 5.1 in Appendix 1 Sec. 5 in [8], or simply by performing a direct check, we obtain 

Proposition 1 Assume that G(Y, t) is twice continuously differentiable in t and that 
sup yt |^G(Y, t)\ + -^G(Y,t) < +oo, which is condition (5.1) in [8, App. 1 Sec. 5j. 
Then, the stochastic process 

G(Y(tA, W),t)- J* (jLG(-, s)](Y(sA, W)) + ^(F(sA, W), s^j ds 

is a martingale with respect to $$ t , the filtration induced by (Y(sA, W))o< s <t- 
For Y standing for some instant of Y(tA, W), define 

My) ■= E w) 

/GMaxPolytopes d _ 1 (y) 

where, recall, MaxPolytopes d _ 1 (F) are the (d— l)-dimensional maximal polytopes of Y (the 
I- segments in the two-dimensional case), whereas </>(•) is a generic bounded and measurable 
functional on (d — l)-dimensional facets in W, that is to say a bounded and measurable 
function on the space of closed (d — l)-dimensional polytopes in W, possibly chopped off 
by the boundary of W, with the standard measurable structure inherited from space of 
closed sets in W. Whereas the so-defined is not bounded, we cannot directly apply 
Dynkin's formula (see Appendix 1, Section 5 in [8] for example) to conclude that the 
stochastic process T,$(Y(tA,W)) — J* * LS0(F(sA, W))ds is a ^-martingale. However, a 
suitable localization argument can be applied (see [17] for the details) to show this: 

Proposition 2 The stochastic process 

^(Y(tA,W))- f [ <j>(f)MdH)ds 

is a martingale with respect to $t t . 



/6Cells(y(sA,W)nH) 
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4 A general variance formula 



The general martingale statements from the previous section admit a convenient special- 
ization to deal with second-order characteristics of iteration infinitely divisible MNW- or 
stationary STIT tessellations. Let us fix through this section a compact convex window 
W C lR d with interior points. From now on we will focus our attention on translation- 
invariant face functionals <fi of (d— l)-dimensional facets, regarded as usual as closed subsets 
of W, of the form 

<f>(f)^Yo\ d ^(f)C(n(f)) (2) 

with n(/) standing for the unit normal to / and ( for a bounded measurable function on 
<S d _i. Recall now the definition of £^(Y(tA, W)), introduce the bar notation E^(F(tA, W)) 
by £ (Y(tA, W)) : = £*(Y(fA, W)) - E£ (Y(tA, W)) and put 

A+{Y{th,W)):= I <f>(f)HdH). (3) 

feCells{Y(tA,W)nH) 

Then we have 

Proposition 3 The two stochastic process 

^(Y(tA,W)) and Til(X(tA,W))- [ A^(Y(sA,W))ds (4) 

Jo 

are both ^s t -martingales. 



Proof of Proposition 4. For some instant Y of Y(tA, W) define 

G{Y,t) :=(^(Y)-E^(Y(tA,W0)) 2 , 
so that G(Y(tA, W),t) = Ej(Y(fA, W)). We use now Proposition 2 to check that 

^G(Y(tA, W),t) = -2[£^(Y(tA, W)) - E^(Y(tA, W))]EA^(Y(tA : W)) (5) 

with given by (3). Put now together (1), Proposition 2 and (5) and use localization 
as in the discussion preceding Proposition 2 in [17] with Gn, N — > oo, chosen so that 
(G N (-,-) A N) V -N = (G(Y,t) A N) V — A r , that \G N (;-)\ < N + 1 and that Gjv(-,t) 
be twice continuously differentiable in t, and with the localizing stopping times Tjy = 
inf{t > : (|G(y(*A,W),t)| V|fG(y(*A,W),i)| V|^G(y(*A,W),*)|) > N}. Proceeding 
as there, we readily conclude that 

Ej(Y(tA, W)) - f I <t> 2 {f)A{dH)ds+ 

J0 J[W] f e Cells(Y(sA,W)nH) 
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2 f\f J2 Hm<t>(Y(sA, W)) - ^(Y(sA, W))]A(dH)- 

JO J[W] f eC en S (Y(sA,W)nH) 

[E*(y(sA, W)) - E^(Y(sA, W^))]E^(y(sA, = 

Zl(Y(tA,W))- [ A t p(Y(sA,W))ds-2 f ^(y(sA, W))S^(y(sA, (6) 

is a 3 r martingale, with A (y (sA, W)) := ^(y(sA, W)) - EA (y(sA, W)). We now take 
advantage of the special form (2) of the face functional to conclude that 

= I VoLj_i(# n W)((n(H))A(dH) = const . 

J[W\ 

This implies A^ = and thus, by Proposition 2 and (6), we can complete the proof. □ 

The so-far established theory is now used to calculate the variance of face functionals as 
given by (2) of iteration infinitely divisible random MNW-tessellations y(tA, W) restricted 
to a compact convex window W C W 1 with VoLj(iy) > 0. 

Theorem 1 For arbitrary diffuse and locally finite measures A on H and <fi as in (2), we 
have 

Var(S,(y(tA,W0))= / C 2 (n(H))[ [ i^^MMli dxdy A { dH). 

J[W] JwnH JwnH Iy \[ x y\) 

Proof of Theorem 1. Recall first (3) and note that it implies 

A^(Y(tA,W))= I <f> 2 (f)HdH) 
J ^ /eCciis(ynH) 

= / C 2 (n(i/)) / / l[x,y are in the same cell of Y n H]dxdyA(dH). 
J[W] JwnH JwnH 

Thus, using (4) and taking expectations of both sides yields immediately 

Var£^(y(tA,W0)= / / ( 2 (n(H))[ [ P(z,yare... 

Jo J[W] JHnwJHnw 

... in the same cell of y(sA, W) H H)dxdyA(dH)ds. (7) 
Taking into account that 

F(x, y are in the same cell of y(sA, W) D H) — exp(— sA([xy])), 
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which follows Theorem 1 [17] together with standard properties of Poisson hyperplane 
tessellation, and using (7), we end up with 

Var(£^(y(*A, W))) = I f ( 2 (n(H)) [ [ exp(-sA([xy]))dxdyA(dH)ds 

Jo J[W] JwnH JwnH 

= / c2 < a < ff »/ / 1 "l ( fM) M)) <to » A(dg) - 

J[w] JwnH JwnH A vF2/JJ 

which completes our argument. □ 

For general hyperplane measures A this cannot be simplified further. However, in the 
special case, where A is the unit-density isometry-invariant measure A iso , tools from inte- 
gral geometry become available to evaluate the integral further. 



5 Exact variance expression for the isotropic STIT 
tessellation 

For the stationary and isotropic case A = Ai SO we want to evaluate the variance expression 
from Theorem 1 further in the special case = VoLj-i, i.e. when £ = 1. To simplify the 
notation we will write from now on Y(t) instead of F(tA iso ). 

Theorem 2 For the stationary and isotropic STIT tessellation Y(t) with surface intensity 
t > we have 

2K d-i t » T I. 
d — 1 f f 1 — e in <i 
Var(Vol d _ 1 (y(*,W))) = -— / / - 2 dxdy (8) 

^ JwJw If 



d(d-l)K d f°° . , H „ / -^±tr x , 

1 Iwi^y 1 - e ^ dr, (9) 



where W is a compact and convex subset ofM d , and where 7vj/( r ) = Js i Vohj(W fl (VF + 
ru))t , d_i((i«) is the isotropized set-covariance function of the window W . 

The key to Theorem 2 is a general integral-geometric transformation formula of Blaschke- 
Petkantschin type, which is interesting in its own right. 

Proposition 4 Let W C M. d be compact and convex and let g : W x W R be a non- 
negative measurable function. Then 

[ [ [ g(x,y)dxdyAUdH)= {d ' 1)Kd - 1 [ [p^\dxdy. (10) 
J[W] JwnH JwnH dK-d Jw Jw IF — y\\ 
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Proof of Proposition 4. First, we use the affine Blaschke-Petkantschin formula [15, 
Thm. 7.2.7] with q — 1 to obtain for any non-negative measurable function h : (R d ) 2 — > M 



h(x,y)dxdy = 



mh(x,y) \\x - y\\ d 1 £ L {dx)£ L {dy)dL, 



where £ is the space of lines in R d with invariant measure dL, i.e. the affine 1-dimensional 
Grassmannian in M. d and ^ is the the Lebesgue measure on L with normalization as 
specified in [15, Thm. 13.2.12]. Taking now 

h(x,y) = l[x E W]l[y E W] \\x - y\\ k g(x, y) 

for some k > —d and another non-negative measurable function g : W x W — > M d we 
obtain 

/ / \\ x ~ V\\ k 9(x,y)dxdy = ^— : f f f \\x - y\\ d ~ 1+k g{x,y)£ L (dx)d£ L {dy)dL. 
Jw Jw 1 Jc JwnL JwnL 

For k — — 1 this yields 

/ / if X,V ,, dxdy = / / / ||rr - y\\ d ~ 2 g(x,y)£ L (dx)£ L (dy)dL. 
Jw Jw \\ x ~~ y\\ 1 J c JwnL JwnL ^12) 

We replace now in (11) for = 0, W by W C\ H for some fixed hyperplane if and d by 
d — 1 and get 

/ / g(x,y)dxdy 
JwnH JwnH 

^ </ VKnffnL JwnunL 

where by £ ff we mean the 1-dimensional affine Grassmannian restricted to H with invariant 
measure dL H (this is the set of lines within hyperplane H). Averaging the last expression 
over all hyperplanes H and using the fact that A iso (dH)®dL H = dL, see [15, Thm. 13.2.12], 
yields 




g(x,y)dxdyA iso (dH) 

'n JwnH JwnH 

(d-l)K d ^ I II I \\ X -yf-Z g ( Xj y)l L ( dx )l L ( d y) dL H Aiso ( d H^ 



H J C H J wnHnL JwnnnL 



= i^htl f f f \\x-y\\ d - 2 g(x,y)£ L (dx)£ L (dy)dL. (13) 
1 J c JwnL JwnL 

By comparing (12) and (13) we finally conclude 

III g(x,y)dxdyA iso (dH)= {d -} )Kd - 1 f f ^\dxdy, 
J[W] JwnH JwnH dKd JwJw\\ x ~y\\ 

completing thereby the proof of the proposition. □ 
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Proof of Theorem 2. In view of the general formula from Theorem 1, we take 



g(x,y) 



1 - exp(-tA iso ([xy])) 1-e 



-t x- 



A iS o(ND ^\\x-y 



where the equality is a simple consequence of the mean projection formula from integral 
geometry, see [15], Thm. 6.2.2 with q = j = d — 1 there. Thus, upon applying the 
transformation formula (10) we conclude the following identity for Var(Sv i d _ 1 (^(t, W))) = 
VartVoWY^WO)): 



f f f d— \ff\ — e iKd ^ x y " 

/ / / g(x,y)dxdyA iso (dH) = — — / / p dxdy 

J [w] J Hnw J Hnw 1 Jw Jw \\x — y\\ 



d(d — \)Kd 



lw( r ) ^ r dr 

/°° / 2K d-i \ 

lw( r ) r ' d ~ 3 ( 1 ~ e ~ tT ) dr, 



2 

where we have passed to d- dimensional spherical coordinates. □ 

In the special case W = B^, the isotropized set-covariance function 7^3 (r) takes the 
form 



^(i-l + i) :0<r<2i? 



:r>2R 
and the variance integral can be evaluated in a closed form: 

4tt 2 

Var(Vol 2 (F(t, B%))) = — {t 2 R 2 {\2 - 8tR + 3t 2 R 2 ) + 24(1 + tR)e~ tR - 24) . 

3t (14) 

The same closed form cannot be obtained for d = 2, since 7^2 (r) has a more complicated 
structure, i.e. 

7 R 2 (r) = 2R 2 arccos ( — ) V 4R 2 — r 2 

R ' \2R) 2 

for r between and 2R and 7^2 (r) = for r > 2R. Unfortunately, the resulting integral 

can in this case not further be simplified. 



6 The variance in the asymptotic regime 

Another important task in our context is to determine for fixed t the large R asymptotics 
of the variance Var(VoLj_i(Y (t), W R )) for the family of growing windows W R — R- W, R — > 
00, with W as in the previous section. Writing ~ for the asymptotic equivalence of func- 
tions, i.e. f(R) ~ g(R) iff f(R)/g(R) — > 1 as R — > 00, we have 
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Theorem 3 For d = 2, 

Var(Voli(y(t,WOi)) ~ 7rVo\ 2 (W)R 2 \ogR, (15) 
whereas for d > 3 we have 

Var(VoWy(*, W R ))) ~ i*^) Izl^w) (i 6 ) 

wift ^(W) fremg the 2-energy ofW, see [9, Chap. 8], given by 

E 2 {W)= [ [ \\x-y\\- 2 dxdy. (17) 
Jw Jw 

In particular, this establishes weak long range dependencies present in stationary and 
isotropic STIT tessellations Y(t). In the planar case, these dependencies are rather weak 
in that 

Vai(Voh(Y(t,W R ))) „ 

vrlogi? ->■ oo, 

Vo\ 2 (W R ) 

as R — > oo. For d > 3 these dependencies are much stronger, as the variance of the total 
surface area grows asymptotically like R 2 ^ 1 ^. 

Proof of Theorem 3 Formula (15) can be established by using (9), the relation <y w ~ 
VohiWn) = R 2 VohiW), valid uniformly for arguments r = 0(R/ log R), the observation 
that ^ Wr -> for r = Q{RlogR), together with the fact that Jq (R) {1 - e~ cr )^ ~ logR, 
c > 0, as soon as log L(R) ~ log-R, and the scaling property of STIT tessellations: 

7r f°° i dr 

Vnr(Voh(Y(t,W R ))) = r 2 Vex(Voh(Y(l,W tR ))) = -j j w Jl- e ~ r ) — 

~ ixr 2 Vol 2 (W tR ) \og(tR) = it Vo\ 2 (W)R 2 (log R + log*) 
~ nR 2 Vo\ 2 (W)\ogR. 

To see (16), use (8) and again the scaling property of STIT tessellations to obtain 

Var(Vo\ d ^(Y(t,W R ))) = R 2{d ' l) Var(Vol d _i (y(itt, W)) 



2 k 



,d- If f l-e 



d-l 



-Rt\\x-y\\ 

= R^t-± I I dxdy 



2 



w Jw \\x - y\\ 2 



R^-V—E^W), R^oo. 
2 

Observe that this does not extend for the separately treated case d = 2 because there the 
integral in (17) diverges. □ 

It is easily seen that E 2 {-) enjoys a superadditivity property 

E 2 {W X U W 2 ) > E 2 (W X ) + E 2 (W 2 ), w 1 nw 2 = 0, 
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which stands in contrast to (15), where the asymptotic expression is linear in VohiW). 
We will now derive an integral-geometric interpretation for this energy functional. Taking 
g(x,y) = 1 and k = — 2 in (11) yields the identity 



E 2 (W) = I I \\x-y\\~ 2 dxdy = ^ I I I \\x - y\\ d ~ z £ L (dx)£ L (dy)dL 

= (d _ Z _ 2) jf W n Lf-^L = (d _ ^ _ -^ (W) (18) 

with I d ^i{W) being the (d — l)-st chord power integral of W in the sense of [15, p. 363]. 
More precisely, 

I d ^(W) = ^ J Vo\i-\w n L)dL. 

Hence, combining (18) with (16) from above ans using the fact that Id-i(-) is homogeneous 
of degree 2(d — 1), we arrive for d > 3 at 

Corollary 1 The asymptotic variance Var(VoLj_i(Y(£, Wr))), R — > oo, is given by 

Vax(VoU-i(Y(t,W R ))) ~ -i-/ d _ 1 (W fl ) = -i-^V^). 

a — 2 a — 2 

In general, / d _i(VF) cannot further be evaluated. But for in the special case W = Bf we 
have by applying [15], Theorem 8.6.6 (with a corrected constant), 

h-i{Bf) = d2 ^^tl 

Kd-l 

and, thus, the 2-energy of the ci-dimensional unit ball Bf equals 



(d-l)(d-2) (d-l)(d-2) V 2 

In the particular case d = 3 we obtain the value E 2 (Bf) = An 2 , which agrees with our 
explicit variance formula (14) from above. 



7 Pair-correlation function 

It is our next goal to establish a closed formula for the pair-correlation function gd(r) of 
the random surface area measure of a STIT tessellation Y(t), which is a commonly used 
tool in spatial statistics and stochastic geometry to describe the second-order structure of 
a random set. It describes the expected surface density of Y(t) at a given distance r from 
a typical point of Y(t), see [15] or [21] for exact definitions. It is well known (cf. [21, p. 
233]) that the pair-correlation function gd(r) and the variance Var(VoLi_i(Y(i, W))) are 
related by the general formula 

poo 

Var(Vol d _ 1 (y(t, W))) = dn d t 2 / i w (r)(g d (r) - l)dr. 

Jo 
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Figure 2: Comparison of the pair-correlation function for a stationary and isotropic STIT 
tessellation Y(l) and a stationary and isotropic Poisson line tessellation with edge length 
density 1 



Thus, from the explicit variance formula in Theorem 2 the following can directly be de- 
duced: 

Corollary 2 The pair- correlation function Qdij-) of the random surface area measure of 
the stationary and isotropic random STIT tessellation Y(t) is given by 



which was independently obtained by Weiss, Ohser and Nagel by entirely different methods 
and is presented in [22]. However, it should be emphasized though that our original 
approach developed above yields information also on higher dimensional cases. For example 
we have for the spatial case d = 3, 



It is interesting to compare the pair-correlation function of Y(t) from Corollary 2 with 
the corresponding function of a stationary and isotropic Poisson hyperplane tessellation 




Especially for d = 2, gaij-) becomes 
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PHT(t) with the same surface intensity t > 0. The latter will be denoted by g d (r). 
Using Slivnyak's theorem for Poisson processes [15, Thm. 3.3.5] one can easily show that 



PHT(t) /\ _ 1 , (d- l)n d - 



dndtr 



Especially for the planar case d = 2, i.e. for the Poisson line tessellation abbreviated by 
PLT(t), we have g^ 1 '^ (r) = 1 + l/(7rtr). A comparison of <?2(?") and g^^ir) is shown in 
Figure 2. 
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